Nonequilibrium dephasing in Coulomb blockade quantum dots 
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We present a theory of zero-bias anomalies and dephasing rates for a Coulomb-blockaded quantum 
dot, driven out of equilibrium by coupling to voltage biased source and drain leads. We interpret 
our results in terms of the statistics of voltage fluctuations in the system. 
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Nonequilibrium phenomena encountered in steady- 
state quantum transport through nanoscale or meso- 
scopic conductors are presently attracting a lot of inter- 
est. Recent experiments on mesoscopic wires have shown 
that the nonequilibrium distribution function evolves 
from a two-step to a broad single-step distribution as 
the effective interaction strength increases HI. These ex- 
periments have led to rich and nontrivial insights into 
dephasing and interaction physics in general. In the pres- 
ence of a finite current flowing through the system, noise 
will be generated, which then leads to decoherence and 
dephasing processes. This phenomenon has been studied 
in some detail for the nonequilibrium Kondo problem, 
both experimentally [2] and theoretically [3], and a zero- 
temperature decoherence rate T{V) ~ V/\n (V/Tx) was 
reported in the regime V ^ Tk- (Here F > is the 
applied bias voltage, T^ the Kondo temperature, and we 
set fi.= e = /cs = l.) Similar dephasing rates were found 
for a nonequilibrium Fermi-edge singularity problem [3], 
and for a spin-fermion model driven out of equilibrium 
[5 . Related questions are also considered for disordered 
interacting quantum wires (Luttinger liquids) [6J , and in 
the context of electronic interferometry [7] . Nonequilib- 
rium dephasing rates can be conveniently inferred from 
the voltage-induced broadening of the zero-bias anomaly 
(ZBA) [5|, i.e. a smearing of the dips in the energy- 
dependent tunneling density of states (TDoS), see also 
Refs. lUITU]. (Note that in general, dephasing rates will 
depend on the quantity under consideration l6].) 

In this work, we present a theoretical study of the per- 
haps most simple and basic scenario where one can study 
such questions, namely for a metallic Coulomb-blockaded 
quantum dot [IT], with a transport voltage V applied 
via attached source and drain electrodes. Note that the 
often-studied single-electron box (SEE) P^ [Hj [Til [T5] 
corresponds to a single-lead situation plus a capacitively 
coupled gate, where one cannot have steady-state current 
flow. Previous theories for the two-lead setup have only 
considered the IV curve (or related transport quantities) 
in certain limits [121 EZ]. As we demonstrate below, the 
system is sufficiently simple to allow for a quantitative 
study of nonequilibrium dephasing. At the same time, it 
displays rich behavior that can be probed experimentally 
with available setups. 

Nonequilibrium effects are captured by the Keldysh 
formulation [^, TH] , where the total system evolves from 



the initial time —to/2 after which the interaction is 
smoothly switched on, to time +to/2 and back; one fi- 
nally takes the limit to -^ cx). By a sequence of standard 
steps [HI HU US] , we obtain a Keldysh functional integral 
representation as 

(X)^ E fT^i^,Qy''^'^''^^^'''-^^^X[Q,ct.], (1) 
wez'' 

where X[Q„,(j)^] may be any observable expressed in 
terms of the charge on the dot (Q<j(i)) and its phase 
{(j>cr{t)), both defined on the upper/lower (tr = ±) branch 
of the standard Keldysh contour [18 . The phase fields 
obey the boundary conditions, X)<j=± '^4'a-i~to/2) — 
2ttW, J2a=±'^^cyi+'^o/'^) = 0, where the integer W is 
summed over [20]. Here, W is the real-time analogue 
of the winding numbers central in establishing Coulomb 
blockade physics in imaginary-time theories [141 |T5]. The 
charging energy, Ec, of the dot enters the theory through 



/•to/2 

5c = E^/ dt[-E,{Q. 

a ''-to/2 



Qgf + Qadt(l)a], (2) 



where the constant Qg defines the electrostatically pre- 
ferred charge configuration. The coupling to the source 
and drain {a = ±) electrodes, biased by aV/2, respec- 
tively, leads to the tunnel action 12, 16 



Sf 



2 ^ 



dtdt' e^^'^-W L„„,{t-t') e^-^-'^*'), (3) 



where, for simplicity, we assume identical tunneling 
conductances gr for both contacts; the generalization 
to asymmetric cases is straightforward. The complex- 
valued functions L^^ji (t) appearing in Eq. ([3| are 






acr'L'{u}) + {a - <t')lj] , (4) 



where the applied voltage V and the temperature T enter 
only via the real part. 
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Our observable of main interest, the energy-dependent 
TDoS i^(e, V, T) (technically, {—tt~^) times the imaginary 



V ^<= 

9t^ 



E. 



c 



9t^ 



FIG. 1: Classical resistor network equivalent of the system. 
For a discussion, see main text. 



part of the Keldysh retarded Green function) affords the 
representation v = v^ + Vh, where v^ (i^h) denotes the 
contribution of electrons (holes) tunneling onto the dot, 

"^-Re /"dTe"^(l-n)(T)(e'(^-(*+^)^*+(*^)). (6) 

Here, vq is the noninteracting DoS, t is an arbitrary ref- 
erence time, n^r) denotes the Fourier transform of the 
double-step distribution function, n(e) = ^{nf{e+V/2) + 
nf{e — V/2)), and nf{e) = (e'^/"^ -I- 1)""'^ is the Fermi dis- 
tribution function. The hole contribution obtains by ex- 
change (1 — n) ^^ n and (j)± ^^ (pzp- We finally note that 
the current flowing through the dot is given by 



I{V) = 



9t 



de[nf{e~V/2)~nf{e + V/2)] 



Ke) 



(7) 



with I ^ Iq — {gT/2)V in the absence of interactions 
{Ec = 0). The phase factors in Eq. (|6l encapsulate the 
effects of particle interactions. We next explore the effect 
of these phase fluctuations on the TDoS for the case qt ^ 
1, but return to the opposite limit, gx ^ 1, at the end 
of the paper. 

In the "open" limit, (/t 3> 1, charge fluctuations are 
large, and the phase representation obtained by integrat- 
ing out the charge flelds in Eq. (fTl) can be restricted to 
a quadratic approximation in the zero winding number 
sector, W = 0. Interaction effects are then relatively 
weak, as in the corresponding SEB problem |14j . and 
the dot distribution function remains close to the dou- 
ble step n{e). In this limit, the "collective variables" (Un- 
essentially describe the classical non-equilibrium steady 
state of a biased RC circuit, see Fig. [T] To explore this 
connection, we expand the phase functional to second 
order in the Keldysh "classical" and "quantum" fields, 
(pc = ((/'+ + 4>^)/2 and 4>q = 0+ — 0_, respectively. 
A "Hubbard-Stratonovich" transformation of the con- 
tribution of 0{(t)'i) then leads to the partition function 
Z ~ /l?(0c, 0g) ^e'^'^**.[(-ca?-2<,,9.)0.+«]^ ^ ^Ijg^e ^g 

identified E^ = (2C)~^ with an inverse classical capac- 
itance. The Hubbard-Stratonovich noise field ^ with 



(^(t))^ = is correlated as {i{t)i{t'))^ = gxL'it - t'), 
and L'{t) is the Fourier transform of the kernel ((sl). Inte- 
gration over (j)g then constrains fluctuations of the volt- 
age U = dt<j>c (relative to the stationary value V/2) 
to solutions of the semiclassical i?C-Langevin equation 
CdtU{t) + R-^U{t) = ^{t), where R = 1/(25t) is the 
parallel resistance of the circuit in Fig. [11 and the fluc- 
tuating noise current /a* = (2Ai)^^/j dt'^{t'), av- 
eraged over a characteristic time window Ai, gives rise 
to the dc noise power S = 2 At var(/) ^I\. At flnite 
temperatures and zero bias, the classical limit of Eq. (|5|, 
L'{t) = 4TS{t), describes a thermal equilibrium situa- 
tion, where the voltage U{t) relaxes to the Boltzmann dis- 
tribution P{U) ~ exp(— Ct/2/2T), and Johns on- Nyquist 
thermal noise, S — AT{gT /2), is recovered. For high volt- 
ages and low temperatures, T <^V , the noise correlator 
([5]) instead asymptotes to L'{t) — V6{t), which impHes 
shot noise, S = 2FIo, with the expected [21! Fano fac- 
tor F = 1/2: upon increasing V, transport through the 
quantum dot undergoes a crossover from an equilibrium 
thermal to a steady-state nonequilibrium shot-noise dom- 
inated regime. 

Employing the TDoS ([6| as a reference observable, we 
next explore the dephasing influence that this noise has 
on the quantum physics of the system. To this end, we 
use (I)c{t2) — (t>c{ti) — J^^ dt U{t) in Eq. ([6| and integrate 
over (pq. As a result, the (electron contribution to the) 
TDoS assumes the form 



= Re /"dTe""(l-n)(T)<Je'^/*"'^*^=W^ , (8) 



t^e(e) 



where Ue is the solution ^22, to a variant of the previous 
Langevin equation, 

{Cdt + R-^)U, = e + Sc [5{t -t) + 5{t - {t + t))] . (9) 

Equation ((8| affords an intuitive interpretation: the 
TDoS Ve probes tunneling processes where an electron 
enters the dot at time t and leaves at time i-\- t. In the 
present structureless environment (a "dot"), the dynam- 
ical phase controlling these processes obtains by integra- 
tion of the time-dependent voltage on the dot. The latter 
is governed by a superposition of (i) a voltage pulse of 
height Ec upon the entry of the external particle, as de- 
scribed by the second term on the r.h.s. of Eq. ([9l j23j. 
and (ii) the background voltage noise ^ on the dot. Av- 
eraging the exponentiated solution of (|9l over the ^ and 
adding the hole contribution, Vh, we obtain 



vo 47rgT Jo ^ 



cos[(e + aV/2)T] 



(10) 



where non-singular contributions of higher order in gj, 
have been neglected. In Eq. ( 10 ), the first term under the 



integral is the temporal Fourier transform of the distribu- 
tion function, the factor (1— e~^*), fi = {RC)~^ accounts 



for the relaxation of the initial voltage pulse, and the last 
term defines the average noise-action (cf. Eq. dsl 



Sir) 



n^ 



Airgj 



dw 



1 — cos{llJt) 



L'H. (11) 



Note that in this weak Coulomb blockade limit, the TDoS 
is independent of the reference charge Qg. In what fol- 
lows, we consider the limit T = where the infiuence of 
nonequilibrium dephasing on the TDoS is strongest. 

To logarithmic accuracy, the equilibrium (V — 0) 
limit of Eq. (Ill can be approximated by 5*0 (t) ~ 
{2Tr gT)~^ln{UTy The resulting TDoS is symmetric in 
e and displays a characteristic dip at e = - the ZBA - 
which vanishes in the limit Ec ^ 0, and also for gj^ — > oo. 
At finite voltage, the TDoS remains symmetric, while the 
double-step profile of the distribution n(e) implies a split- 
ting of the e = ZBA into two minima at e = ±^/2, see 
Fig. |2] The strict positivity of S{t) - S'o(t) > 0, pre- 
viously identified as a manifestation of enhanced noise 
levels for V ^ 0, then causes a suppression of the ZBA. 
At the same time, its line shape broadens. In the limit 
V — > cx), the linear growth 5(t) ~ Vt implies a van- 
ishing ZBA. Similar phenomena were recently studied [8 
for tunneling into a 2D diffusive metal. Eq. (10) may 



serve to define a nonequilibrium dephasing rate T{V) 
from the voltage-induced broadening of the ZBA dips at 
e = ±V/2 0. Writing Se = e - V/2 with |<5e| < n, the 
ZBA dip is described by Si^{Se) = ^q — i^(e). Parameteriz- 
ing small deviations off the dip at V/2 in terms of the de- 
phasing rate [24,, [Sv{0) - Siy{5e)]/Si^{0) ~ ^[5e/T{V)]^, 
Eq. ([To| yields 



T^iV) = 



/-^e-^M(l-e-"^) 
J^dTTe-si-){l-e-^^^)' 



(12) 



Analytical results for T{V) can then be extracted from 
Eq. (12 1 in various limiting regimes, while the full curve 
is obtained numerically, see Fig. [2] 

Let us first discuss the asymptotic regime of weak in- 
teractions, Ec -^ 0. (Technically, this means that we 
take the limit fi ^ prior to Se -^ 0.) Eq. (12 1 is 
then dominated by contributions ilr <ti 1, which im- 
phes S{t) w VEcT^/4: and r(y) = y/EcV/2. The gr- 
independence of F refiects the diverging _RC-time: par- 
ticles tunneling onto the dot do have time to realize 
the dissipative nature of their environment. Turning 
to the complementary regime, (Se — > at fixed fi, we 
note that the leading contribution to an expansion of the 
TDoS in g^^ <C 1, obtained by setting 5* = 0, is given 
by the logarithmically singular expression Sv{de)/i'o = 
(87rgr)"^ln[l -I- {fl/6e)^]. Using Eq. ^, this recovers 
the well-known result [THJ [17] for the differential con- 
ductance, {2/gT)dI/dV = 1 - (47r5T)"^ln[l + i^/V)^]- 
However, in order to define the dephasing rate, one needs 
to retain a finite noise action S. Specifically, for V > il, 
S{t) ~ Vr/iAirgT), implying 



V 



1 



47rgT 



ln(l + Angrn/V) 
l-{l + AngTn/V)- 



1/2 



(13) 
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FIG. 2: (Color online) Dephasing rate vs voltage (at T = 0) 
in units of the inverse RC time, Q, — igrEc, for different 
values of qt ■ Curves were obtained by numerical integration 
of Eq. ( 12 1. Dotted (al) and dashed (a2) curves represent the 
analytical predictions in Eqs. (131 and (141, respectively; for 
(a2), c = 0.25 in Eq. ^. The insets show the TDoS iy{e)/iyo 
at r = as a function of e for several V. Left inset: from 
numerical integration of Eq. ( 10 1 for gr = 5. Right inset: 
from Eq. (fTel) for gr < 1 and Qg = 0. 



while for fie" '"^^ <V <n f25], we estimate 



V 



V^T^gr 



(v/n) 



-l/(47rgT) 



(14) 



where c is a numerical prefactor of order unity. Impor- 
tantly, in all parameter regimes, we find T(V) < V, and 
the double peak structure in the nonequilibrium TDoS 
is reasonably well resolved. As a function of V, the 
relative strength of the dephasing rate (i.e. the ratio 



F/y) increases with decreasing V. Furthermore, Eq. ( 13 1 
shows that for very large voltage, V > gr^, the rate 
T(V) ^ V/gT is directly determined by the shot noise 
5 = 
and 



2FIo discussed above. The analytical results (13 1 
(ITll) describe the numerical solution to Eq. ( 12 1 



rather well, as indicated in Fig. [2j With increasing gr, 
the dephasing rates are gradually suppressed, and for 
F = 0, as expected, the dephasing rate is zero. 

Finally, we briefiy turn to the opposite case of strong 
Coulomb blockade, gr <^ 1, where the dual charge repre- 
sentation is appropriate. Expanding e*'^'"" in Eq. (IT|) into 
a Taylor series, we can integrate out the phase fields (j)a-. 
Now the e^*'^"'^*-' factors appearing in the tunnel action 
^ can be interpreted as charge raising or lowering oper- 
ators, i.e. they change the corresponding Qait) variable 
by ±1. In order m of the tunnel expansion, we thus have 
2m jumps at times tj, t^- (j = 1, . . . , to), with jump direc- 
tions cTjiCr' = ±1. The charge fields are then expressed 
in terms of these jump times and directions, Qa{t) = N — 

^Y:T=i [s<ra,eit-tk)~S„„>^eit~t'^)\ , where e(t) is 
the Heaviside function and the ly-summation in Eq. (fTl) 
implies the boundary condition Q±{—to/2) = N E Z. 



The charge representation of Eq. (fTl) is then given by [5B] 

^-E^i^) E /H'^Mt; (15) 



m— 






fe=i 



Away from charge degeneracies, i.e. for Qg not half- 
integer, Eq. ( 15 1 admits a solution in terms of a 



■'noninteracting-blip approximation" (NIB A) |27j . where 
Z is dominated by short "blips" of length \tj — t'A 
during which Q+(t) ^ Q^(t), with long time inter- 
vals between subsequent blips where Q+{t) = Q-(t). 
Defining the optimal charge state Q G Z according to 
Qg-l/2 <Q <Qg + 1/2, the NIBA transition rate from 
Q± — Q to any other charge state vanishes for T — Q 
and F < 2£;c|Q - Qg + 1/2|. In this regime, the TDoS is 
thus given by (see Fig. [2]) 



K£) 



E 



s — ±,Oi. — ±: 



'-Q(e+[s + 2{Q-Qg)]E,+ '^ 



(16) 
which recovers the equilibrium result of Ref. [15 for 
V = Q. Note that v{e) is generally asymmetric; Eq. (l7| 
implies well-known Coulomb blockade expressions for the 
IV curve. Vanishingly small (noise) current levels imply 
that nonequilibrium dephasing is strongly suppressed in 
the Coulomb blockade regime and may leave traces only 
beyond NIBA. The near-degeneracy case, on the other 
hand, is related to a many-channel variant of the Kondo 
effect [13], where frequent transitions between two charge 
states occur. NIBA does not hold in this regime, but a 
real-time renormalization-group approach [51 can be con- 
structed and will be described elsewhere. 
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